Lane reduction in driven 2d-colloidal systems through microchannels 
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The transport behavior of a system of gravitationally driven colloidal particles is investigated. The particle 
interactions are determined by the superparamagnetic behavior of the particles. They can thus be arranged in a 
crystalline order by application of an external magnetic field. Therefore the motion of the particles through a 
narrow channel occurs in well-defined lanes. The arrangement of the particles is perturbed by diffusion and the 
motion induced by gravity. Due to these combined influences a density gradient forms along the direction of 
motion of the particles. A reconfiguration of the crystal is observed leading to a reduction of the number of lanes. 
In the course of the lane reduction transition a local melting of the quasi-crystalline phase to a disordered phase 
and a subsequent crystallization along the motion of the particles is observed. This transition is characterized 
experimentally and using Brownian dynamics (BD) simulations. 



Pedestrians on a walkway or in a pedestrian zone generally 
move in a well-ordered fashion once the number of people 
per area exceeds a certain threshold The reason for this 
behavior is that people moving in one direction try to avoid 
crossing paths with people moving in the opposite direction. 
This can be most easily achieved, if so-called lanes are formed 
in the motion of the individuals. This phenomenon has been 
observed theoretically in a large number of systems Ql 01 ■ 
Under certain circumstances it can be shown that lane forma- 
tion is also favored in many-particle systems or 
even in the motion of animals |i8J . Most reports of lane for- 
mation have been based on simulations or theoretical calcu- 
lations. Only recently lane formation could be demonstrated 
experimentally in a three dimensional system of oppositely 
charged colloids driven by an external electric field |9|] . 

In this work we report on studies of the transport behavior 
of colloids in a quasi two-dimensional (2d) setup. The col- 
loids are superparamagnetic, therefore the interaction energy 
can be continuously tuned by the application of an external 
magnetic field. The particles are driven through a narrow con- 
striction (channel). In addition to the analogies mentioned ear- 
lier, such a system resembles the classical case of systems like 
a quantum point contact in mesoscopic electronics 
These contacts show transport in electronic channels due to 
quantization effects. A classical version of a similar scenario 
can be built on a liquid helium surface, which is loaded with 
charges. In this system the formation of channels has been 
reported as well 11211 . The main advantage of the use of su- 
perparamagnetic colloids instead of electrons is given by the 
size of the colloids, which can be easily monitored in a stan- 
dard video microscopy setup. All relevant parameters can be 
gathered from the configuration data. 

We compare the experimental results gained from video mi- 
croscopy with BD simulations of particle flow through con- 
strictions under very similar conditions. The results of both 
cases show qualitative and quantitative agreement. The main 
focus of this discussion will be on the lane reduction transi- 
tion, giving possible scenarios for the origin of this transition 
and describing the behavior of the particles in the transition 
region and its close vicinity. 

Two particle reservoirs and a connecting channel are de- 



fined on a lower substrate using UV-lithography. SEM pic- 
tures of the channel setup and of the channel entrance together 
with some dried particles inside and outside of the channel 
can be seen in Fig. Q The channel (60 itm wide and 2 mm 
long) is filled with a suspension of superparamagnetic parti- 
cles, which are commercially available (Dynal, particle diam- 
eters a = 4.55 |J.m, suspended in water). Gravity confines 
the colloidal particles to the surface of the channel due to the 
density mismatch between the particles and the liquid. So the 
system is quasi two-dimensional, as long as the magnetic in- 
teractions do not lead to out-of-plane motion of the particles. 
Thus for an applied uniform magnetic field perpendicular to 
the monolayer the particle interaction is purely repulsive, and 
its strength at distance is given by 



(1) 



with the magnetic dipole moments M = XesB of the parti- 
cles, which are proportional to the external magnetic field B. 
If no boundary conditions are imposed on the system, 2d 
phase transitions from a liquid phase to the hexatic phase 
and finally a crystalline phase can be observed as a func- 
tion of the particle interactions fl3ll . with boundary conditions 
the behavior becomes more complex In transport 

through narrow channels the system orders in lanes at inter- 
mediate strength of the interactions. The importance of the 
pair-interaction can be characterized by the dimensionless in- 
teraction strength fisll 



(2) 



where p denotes the number density of the particles, ks the 
Boltzmann constant, and T the temperature. 

Even at moderate tilts of the experimental setup the par- 
ticles cannot escape the channel and gravitationally induced 
transport between the two reservoirs has to take place in- 
side the channel. At this moderate tilt the velocity of non- 
interacting particles depends linearly on the tilt. A small 
barrier surrounding the channel prevents additional particles 
from falling into the channel. Thus the number of particles 
is conserved in the system reservoirs-channel. The video 
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FIG. 1: SEM-pictures of the channel: Shown are an overview of the 
channel and a zoom to the region of the channel entrance. Some 
dried particles inside and outside of the channel can be seen as well. 

microscopy setup allows for moderate tilting without losing 
optical alignment. Before starting the experiments the sys- 
tem is set up completely horizontal. We use configurations 
where the particles are either all confined in one reservoir 
or equally distributed along the channel and in both reser- 
voirs. The susceptibility of the colloidal particles is Xcff = 
3.07- 10"" Am^/T. Thus an external magnetic field B = 
0.24 mT corresponds to T w 2.5. 

The Brownian dynamics (BD) simulations are based on an 
overdamped Langevin equation. This approach neglects hy- 
drodynamic interactions as well as the short-time momentum 
relaxation of the particles. Both approximations are fully jus- 
tified in the current experimental context. Typical momen- 
tum relaxation times are on the order of 100 \\s and therefore 
much shorter than the repetition rate of the video microscopy 
setup (10 s). The colloidal trajectories r,i(t) = {xi{t),yi{t)) 
[i = 1, . . . , N) are approximated by the stochastic position 
Langevin equations with the friction constant ^ 

= -V, y: y^An,) + Fr + F.(t). (3) 

The right hand side includes the sum of all forces acting 
on each particle, namely the particle interaction, the con- 
stant driving force F™* = mgsm{a)x for the inclination 
a and the random forces Fi{t). The latter describe the col- 
lisions of the solvent molecules with the ith colloidal parti- 
cle and in the simulation are given by random numbers with 
zero mean, (Fi(t)) = 0, and variance {Fia{t)Fii3{Q)) = 
2kBT£^5{t)Sij5af3- The subscripts a and /3 denote the Carte- 
sian components. These position Langevin equations are inte- 
grated forward in time in a Brownian dynamics simulation us- 
ing a finite time At and the technique of Ermak 11161111711 . Par- 
ticles are confined to the channel by ideal elastic hard walls in 
y-direction. The channel end is realized as an open boundary. 
To keep the overall number density in the channel fixed, every 
time a particle leaves the end of the channel a new particle 
is inserted at a random position (avoiding particle overlaps) 
within the first 10% of the channel, acting as a reservoir 




FIG. 2: (a) Video microscopy snapshot of colloidal particles mov- 
ing along the lithographically defined channel, (b) Simulation snap- 
shots for a channel (692 x 60 (xm, F ~ 2.5) with ideal hard walls 
(573.3 X 45 \im, F — 115), (c) the same as in (b) with the particles 
at the walls (marked green) kept fixed (573.3 x 45 (xm, F = 902). 
The blue rectangles mark the lane transition region. 

Starting from a random particle distribution within the 
channel, we first calculate an equilibrium configuration 
^poxt _ 0) of a closed channel with ideal hard walls. Af- 
terwards we apply the external driving force and allow the 
system to reorganize for 10^ time steps, before we evaluate 
the configurations. The time step At = 7.5 • 10~^tB is used, 
with ts = ^a^ /ksT being the time necessary for a parti- 
cle in equilibrium to diffuse its own diameter a. We choose 
^ = 'iirrja, with jy denoting the shear viscosity of the water 
The simulations are done with 2000 — 4500 particles, for a 
channel geometry of L^. — SOOcr and Ly = (9 — 12)(t, and 
r-values ofr « 70- 950. 

A typical snapshot from the experiment of the particles 
moving along the channel is shown in Fig.|2ja). Similar snap- 
shots we get from simulations with co-moving (Fig.|2jb)) and 
fixed boundary particles (Fig.|3c)), i.e., the velocity is kept 
to zero for the particles at the channel wall. In most regions 
of the channel the particles are placed in a quasi-crystalline 
order. This behavior is due to the strength of the particle inter- 
actions caused by the external magnetic field (high F-values), 
which leads to quasi-crystalline behavior in unbounded sys- 
tems as well. The formation of this order naturally gives rise 
to the formation of lanes in the motion of the particles along 
the channel. A similar layering phenomenon has been ob- 
served in channels under equilibrium condition |3i]. Addition- 
ally to this lane formation we observe, both in experiment and 
in simulation, a decrease of the number of lanes in the direc- 
tion of motion. In between both regions therefore a region 
exists in which the particles cannot be well-ordered. This re- 
gion is called the lane-reduction zone. In Fig.|2]these regions 
have been marked. 

The reduction of the number of lanes originates from a den- 
sity gradient along the channel. The local particle density in- 
side the channel is shown in Fig. |3ja) and (b) together with 
the particle separations in x- and y-directions. In the exper- 
iment (Fig. |3ja)) the density decreases monotonically along 
the direction of the motion of the particles by about 20%. The 
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FIG. 3: Local lattice constants and dy and local particle density (a) in the experiment and (b) in the BD simulation. The results are obtained 
for the systems of Fig.|2{a) and (b) respectively, (c) Potential energies per particle of different lane configurations as a function of the particle 
density. The dots mark the perfect triangular lattices for 5, 6 and 7 lanes. Also shown are parts of the configurations with 7 and 6 lanes at the 
intersection point, (d) Plots of the lane order parameter for the configuration snapshot of Fig.|2jb). 



average density in the channel shows fluctuations on the order 
of 10% as a function of time. The total increase in density, 
however, is less than 3% during the total time of the experi- 
ment. We therefore argue that the density gradient is formed 
in a quasi-static situation. This argument is confirmed by re- 
sults of BD simulations (Fig. |3b)), where the corresponding 
decrease of the particle density is observed. 

The particle separations of neighboring particles in x- and 
y-direction are used to calculate the local lattice constant d. 
Due to the density gradient along the channel, the crystal is 
not in its equilibrium configuration at all points along the 
channel. Thus the local lattice constant d^, calculated from 
the particle separations in x-direction, can deviate from the 
local lattice constant dy, calculated from the particle separa- 
tions in y-direction. At the left end of the channel, dx is larger 
than dy, indicating that the crystal is stretched along the x- 
axis. At the point of the lane-reduction transition the crys- 
tal changes back to a situation, where dx is smaller than dy. 
This is achieved by decreasing dx and increasing dy by about 
20% simultaneously. Due to this behavior no non-monotonic 
change in the local density can be observed at the position 
of the lane-reduction transition. The behavior of the system 
shows that the stretching of the crystal before the transition 
causes an instability towards decreasing the number of lanes. 
This decrease compresses the crystal along the x-direction, 
but apparently lowers the total energy of the system. 



This can be qualitatively confirmed by the following rough 
estimation: Starting from an ideal triangular configuration 
with a given number of lanes (ni) in a channel of fixed width, 
we calculated the potential energy per particle for different 
particle densities by just scaling the channel length. Plots of 
these energies per particle for different values of ni as function 
of the particle density are shown in Fig.|3lc). They show clear 
intersection points, indicating that for a stretched configura- 
tion with ni lanes in x-direction it can become energetically 
more favorable to switch to a compressed configuration with 
{ni — 1) lanes. Also equilibrium BD simulations (F^"^' = 0) 
of closed channels with non-parallel walls result in a density 
gradient with decreasing channel width and show lane transi- 
tions. A snapshot is shown in Fig.|3c). 

The region of lane-reduction can be well localized by an ap- 
propriate local order parameter We therefore divide the chan- 
nel of width Ly into several bins in .x-direction each contain- 
ing ribin particles and evaluate for different number of lanes 
(ni) the lane order parameter 



lane,ni 



1 "bi" . 2^{„ 



- Vj 



(4) 



which is unity for ni particles distributed equidistantly across 
the channel width starting at y = 0. As can be seen in 
Fig. |3ld), the lane order parameter exhibits a clear disconti- 
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FIG. 4: Delaunay triangulation of the particles moving in the chan- 
nel. The particles are coded according to the number of nearest 
neighbors they have. Open circles mark the bulk particles with 6 
nearest neighbors and the edge particles, symbol x corresponds to 
fivefold symmetry, and symbol V to sevenfold symmetry, (a) Ex- 
periment: In order to minimize the effects of fluctuations on a short 
time scale, 50 images have been averaged, (b) BD simulation for 
a channel with parallel walls, (c) Equilibrium BD simulation for a 
channel with non-parallel walls (ALy = la between both channel 
ends, r = 28.0). 



the bulk lanes have to change the lane, causing the transition 
to move a little bit in direction of the flow. A particle chang- 
ing into the edge lane can neutralize the defect of the transition 
locally. This causes a reconfiguration of the crystal, which in 
turn gives rise to repositioning of the lane reduction zone back 
to a region of higher density, i.e., back to the closest density 
mismatch between bulk and edge particles. 

To summarize we have shown the formations of lanes and 
the occurrence of a lane-reduction transition in a 2d system of 
superparamagnetic colloids both experimentally and by Brow- 
nian dynamics simulations. The lane formation is induced by 
the repulsive particle interactions and the confining potential 
of the channel. Due to those two factors, a crystal forms in the 
channel. This crystal is stretched along the direction of the 
motion of the particles as a response to the boundary condi- 
tions at both ends of the channel and the particle interactions. 
Due to this stretching a lane-reduction transition occurs in the 
channel, at which the crystal locally disorders. After the tran- 
sition the crystal is compressed in the direction of motion. 

We gratefully acknowledge the support of the SFB 513, the 
SFB TR6 and the NIC, HLRS, and SSC. 



nuity at the position of the lane-reduction. The width of the 
transition region is very sharp and usually occurs within 2 — 3 
lattice constants. The local orientational order parameter ^fg, 
which is often used for 2d systems HjIi . is not so significant 
for this system, as it is very sensitive to any perturbation of 
the sixfold symmetry. In addition to the energy instability the 
transition requires a seed to occur, as will be demonstrated in 
the following discussion. 

Important information on the nature of the lane-reduction 
transition can be governed by the so-called Delaunay trian- 
gulation and the results are shown in Fig. |3 It reveals that 
the system is perfectly crystalline left and right of the transi- 
tion region. The transition is marked by a single defect only. 
In addition to this, periodic defects can be seen close to the 
walls. As can be observed in Fig. |2] already, the particle den- 
sity in the lanes at the walls (edge particles) is higher than for 
"bulk" particles. Due to the long range nature of the repul- 
sive particle interactions it is energetically more favorable to 
place additional particles in the lane at the wall than in a bulk 
lane t^(7(]. This density mismatch between edge and bulk par- 
ticles causes periodic defects, which are confined to the walls. 
The lane reduction is triggered by such a defect in addition to 
the slowly varying particle density along the channel, which 
makes the crystalline system unstable against perturbations. 
The defect causes a local phase transition from crystalline to 
disordered behavior 

Since the position of the lane reduction zone is mainly de- 
termined by the density gradient, its location remains stable 
with time on average. A more detailed analysis reveals, how- 
ever, that the transition zone oscillates back and forth around 
this average position. At the transition the driven particles in 
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